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m (54) Tide: EFFICIENT IMAGE RECONSTRUCTION ALGORITHM FOR VARIABLE PITCH SPIRAL COMPUTED TOMOG- 
^ RAPHY 

00 

O (57) Abstract: Reconstructing images of objects spirally scanned with two-dimensional detectors with a novel algorithm under 

^-j. a variable pitch(nonconstant speed), where the object is not restricted to moving at a constant velocity. The object can move at 

O variable speeds (increasing, decreasing, combinations thereof) during the scan of the object. The image reconstruction process is 

O proven to create an exact image of the object under the ideal circumstances. The algorithm can have a convolution-based FBP 

^ (Filtered Back Projection) structure and works very efficiently. The algorithm uses less computer power and combines the benefits 
of Exact Algorithms and Approximate algorithms. An object can be moved at a nonconstant speed through a rotating source and 

^ oppositely located detector. Additionally, at least one source and oppositely located detector can be mounted on a coil stand for 

j^- generating the spiral « 
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. EFFIGIES ALGOMm 
: VARIABLE prfCH SPIRAL COMiPUTED TOMOGRAPHY 

. . .This invention is a Contmuation-In-Part of United States Patent Application 
Serial No. 10/728,136, fildd December 4, 2003 vviich claims the benefit of priority to* 

\ : . O)nfmiwtion-In-Pit of United Stat^ Patent Application- Serial No. 10/389,534 filed 
. Mardi 14, 2003 whicli is S Cb^^ io/389,090 filedMarcb 

Continuatibn-lxi-Part of Serial ko. Wl43,160 filed May 10, 2002 now U.S. Patent 
6,574,25*9, \^ch"claims the benefit of pnori^to U.S. Provisional Application 
ip 60/31^827 filedAugust 16, 2001; • ** ' 

FiELD OF INVENTION 
. \ This invention relates to computer tomography, and in particular to processes 
and systems for reconstructing three-dimensional images from the data obtained by a 
1 5 variable pitch spiral scan of an object, such as When the object moves at a variable 
speed, while the x-ray source rotates around the object. 

BACKGROUND AND PRIOR ART 
Over the last thirty years, computer tomography (CT) has gone from image 
20 reconstruction based on scanning in a slice-by-slice process to spiral scanning. From 
the 1970s to 1980s the slice-by-slioe scanning was used. In this mode the incremental 
motions of the patient on die table through the gantry and the gantry rotations were 
performed 6neaito another; ' 
. rotations,.the trajectory of the x-ray source around die patient was circular. Pre- 
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'■ . : selected slices. thtoughlhe pifienttiave been reconnected u$ing the dpfia obtained by • 
such circular scans,' * ! 

From the mid 1980s to present day, spiral type scanning has become the 
' preferred process for data collection in CT. Under spiral scanning ia table with the 
c .S^ .-^atfent coptinuQUsly mov.es. at ft pqnstapt speed #)ipug|i th[e gantry thtf is,cpntinRO,usly 
• Rotating ab^irt the table. : At ;fu^.^ detectors!, 
which receive data in<)ne.dimension (ft sjagterovv bf detectors).JLater, two- 
dimensional detectors, where multiple rpws (two or more rows) of detectors sit next to 
one another, have been ihtii>duced r ^In CT.theiS have been significant problems for 
."10-. image' reconstruction especially forWb^timeMional detectors.; In what follows the 
data provided by the two-dimensional detectors will be referred to as cone-beam (CB) 
. -data or CB projections. 

Fig. 1 shoWs a typical arrangement of a patient on a table that moves at a 
, jconstant.speed within a rotating gantry having an x-rfty tube source and a detector 
15 array, where cone beam projection data sets are received by the x-ray detector, and an 
image reconstruction process takes place in a computer with a display for the 
reconstructed image. 

For three-dimensional (also known as volumetric) image reconstruction from 
the data provided by a spiral scan with two-dimensional detectors, there are two 
20 known groups of algorithms: Exact algorithms and Approximate algorithms, that each 
have known problems. Under ideal circumstances, exact algorithms can provide a 
replication of an exact image. Thus, one should expect that exact algorithms would 
"•#-.*. -y-^jyjj^g jjfijjg^ 6f g^'^^Hy W^mder non-id&l (fhafc is,:realistic). circumstances. 
Hovyever, exact algorithms can be known to take many hours to provide an image 
25 reconstruction, and can take up great amounts of computer power when being used. 
These algorithms can require keeping considerable amounts of gone beam projections 
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in memory; Addition;aUy, : *s6mie exact ajgbridiins^ i^uire large defector arrays to ' • 
be operable and cah.have limits on 1hfe size of the patient being scanned. 

Approximate algorithms possess a filtered back projection (FBP) structure, so 
* they can produce an image yeiy efifilcientjy.and using less computing powgr than Ejcact 
<- -.5 : ; -..algbrifta]^. eye^urider .tlie jjdeal circumstances they produce. an . ...... 

. : ^ Iri 
particular,; Approximate ^gjorithhls caii; create artifacts, which; aje.false features in an. . 
: . image. Under certain circumstances th^ artifacts could be quite severe. 

■ : To date, there are no known ^oritftnts that can combine the beneficial 
> 10 attributes of Exact and Approximate algorithms into a single algorithm that 

. of replicating an exact image under the ideal circumstances, uses small amounts of 
. • computer power,.and reconstructs the exactimages in an efficient manner (i.e., using 
the PBP'Structure) in the case of variable pitch spiral scanning. Here and everywhere 
. below by the phrase that the algorithm of the invention reconstructs an exact image we 
15 will mean that in theory the algorithm is capable of reconstructing an exact image. 
Since in real life any data contains.noise and other imperfections, no algorithm is 
capable of reconstructing an exact image. 

Image reconstruction has been proposed in many U.S. Patents. See for 
example, U.S. Patents: 5,663,995. and 5,706,325 and 5,784,481 and 6,014,419 to Hu; 
20 5,881,123 and 5,926,521 and 6,130,930 and 6,233,303 toTam; 5,960,055 to 

Samaresekera et al.; 5,995,580 to Schaller; 6,009,142 to Sauer; 6,072,851 to Sivers; 
6,173,032 to Besson; 6,198,789 to Dafiii; 6,215,841 and 6,266,388 to Hsieh. 
* v H6wev«r, lioiie of &e^>&ttts bV^bine'aiH of the d^clencies to image reconstruction 
referenced above. 

25 
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SUMMARY OF THE INVENTION 
A primary objective of the invention is to provide an improved process and 
system for reconstructing images of objects that have been scanned in a spiral fashion 
with variable pitch(at a nonconstant speed) and with two-dimensional detectors. 
5 A secondary objective of the invention is to provide an improved process and 

system for reconstructing images of objects spirally scanned with variable pitch(at a 
nonconstant speed) that is known to theoretically be able to reconstruct an exact image 
and not an approximate image. 

A third objective of the invention is to provide an improved process and 
1 0 system for reconstructing images of objects spirally scanned with variable pitch(at a 
nonconstant speed) that creates an exact image in ail efficient manner using a filtered 
back projection (FBP) structure. 

A fourth objective of the invention is to provide an improved process and 
system for reconstructing images of objects spirally scanned with variable pitch(at a 
15 nonconstant speed) that creates an exact image with minimal computer power. 

A fifth objective of the invention is to provide an improved process and 
system for reconstructing images of objects spirally scanned with variable pitch(at a 
nonconstant speed) that creates an exact image with an FBP structure. 

A sixth objective of the invention is to provide an improved process and 
20 system for reconstructing images of objects spirally scanned with variable pitch(at a 
nonconstant speed) that is CB projection driven allowing for the algorithm to work 
simultaneously with the CB data acquisition. 

A seventh objective of the invention is to provide an improved process and 
system for reconstructing images of objects spirally scanned with variable pitch(at a 
25 nonconstant speed) that does not require storing numerous CB projections in 
computer memory. 
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An eighth objective of the invention is to provide an improved process and 
system for reconstructing images of objects spirally scanned with variable pitch(at a 
nonconstant speed) that allows for almost real time imaging to occur where images 
are displayed as soon as a slice measurement is completed. 
5 A preferred embodiment of the invention uses a six overall step process for 

reconstructing the image of an object under a spiral scan. In a first step a current CB 
projection is measured. Next, a family of lines is identified on a detector according to 
a novel algorithm. Next, a computation of derivatives between neighboring 
projections occurs and is followed by a convolution of the derivatives with a filter 
10 along lines from the selected family of lines. Next, using the filtered data, the image 
is updated by performing back projection. Finally, the preceding steps are repeated 
for each CB projection until an entire object has been scanned. This embodiment 
works with keeping several (approximately 2-4) CB projections in memory at a time 
and uses one family of lines. 
1 5 Unlike the prior art, the invention is not limited to moving an object at a 

constant speed through a spiral scan. The object can be moved at a nonconstant speed 
through the gantry. 

Other embodiments allow for the object to remain stationary within a spiral 
coil type stand having multiple x-ray sources and oppositely located detectors 
20 arranged along the coil stand which are activated sequentially from different locations 
on the coil stand. Still furthermore, the entire coil stand with fixed plural x-ray 
sources and oppositely located detectors rotates all about the object. 

Still furthermore, the spiral coil stand can contain a single x-ray source and 
oppositely located detector which moves along a spiral track about the fixed object at 
25 constant and nonconstant speeds. Still furthermore, the spiral stand can include coil 
links that are not evenly spaced from one another so that the single x-ray source and 
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opposite located detector pass along the length of the object at different speeds. Thus, 
closely located links allow the single source and detector to pass at a slower rate over 
an object than distantly spaced apart coil links. 

Further objects and advantages of this invention will be apparent from the 
5 following detailed description of die presently preferred embodiments, which is 
illustrated schematically in the accompanying drawings. 



BRIEF DESCRIPTION OF THE FIGURES 

Fig. 1 shows a typical arrangement of a patient on a table that moves within a rotating 
10 gantry having an x-ray tube source and a detector array, where cone beam projection 

data sets are received by the x-ray detector, and an image reconstruction process takes 

place in a computer with a display for the reconstructed image. 

Fig. 2 shows an overview of the basic process steps of the invention. 

Fig. 3 shows mathematical notations of the spiral scan about the object being scanned. 
15 Fig. 4 illustrates a PI segment of an individual image reconstruction point. 

Fig. 5 illustrates a stereographic projection from the current source position on to the 

detector plane used in the algorithm for the invention. 

Fig. 6 illustrates various lines and curves, such as boundaries, on the detector plane. 
Fig. 7 illustrates a family of lines used in the algorithm of the invention. 
20 Fig. 8 is a four substep flow chart for identifying the set of lines, which corresponds to 
step 20 of Fig. 2. 

Fig. 9 is a seven substep flow chart for preparation for filtering, which corresponds to 
step 30 of Fig. 2. 

Fig. 10 is a seven substep flow chart for filtering, which corresponds to step 40 of Fig. 
25 2. 
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Fig. 1 1 is an eigjht substep flow chart for backprojection, which corresponds to step 50 
of Fig. 2. 

Fig. 12 shows an arrangement of scanning an object with a spiral coil x-ray source 
where the object being scanned remains stationary inside. 

5 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Before explaining the disclosed embodiments of the present invention in detail 
it is to be understood that the invention is not limited in its application to the details of 
the particular arrangements shown since the invention is capable of other 
10 embodiments. Also, the terminology used herein is for the purpose of description and 
not of limitation. 



Fig. 1 shows a typical arrangement of a patient on a table that moves within a 
rotating gantry having an x-ray tube source and a detector array, where CB projections 

1 5 are received by the x-ray detector, and an image reconstruction process takes place in 
a computer 4 with a display 6 for displaying the reconstructed image. For the subject 
invention, the detector array is a two-dimensional detector array. For example, the 
array can include two, three or more rows of plural detectors in each row. If three 
rows are used with each row having ten detectors, then one CB projection set would 

20 be thirty individual x-ray detections. 

Fig. 2 shows an overview of the basic process steps of the invention that occur 
during the image reconstruction process occurring in the computer 4 using a first 
embodiment. 

The first embodiment works with keeping several (approximately 2-4) CB 
25 projections in computer memory at a time and uses one family of lines. 
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In the first step 10, a current CB projection set is taken. The next step 20 
identifies a set of lines on a virtual x-ray detector array according to the novel 
algorithm, which will be explained later in greater detail. In the given description of 
the algorithm it is assumed that the detector array is flat, so the selected line can be a 
5 straight tilted line across the array. 

The next step 30 is the preparation for the filtering step, which includes 
computations of the necessary derivative of fee CB projection data for the selected 
lines. 

Hie next step 40 is the convolution of the computed derivative (the processed 
10 CB data) with a filter along lines from the selected family of lines. This step can also 
be described as shift-invariant filtering of the derivative of the CB projection data. In 
the next step 50, the image of the object being computed is updated by performing 
back projection. 

The invention will now be described in more detail by first describing the main 

1 5 inversion formula followed by the novel algorithm. 

Unlike the prior art, the invention can be used with objects that move at 
variable speeds through a rotating gantry. The object can accelerate, decelerate or 
combinations thereof. A slower speed through the rotating gantry can provide 
enhanced images of particular portions of an object as desired. 

20 Experimentation in which the speed of the moving table through a rotating 

gantry was ramped up approximately 25%(twenty five percent) over the course of 
three(3) gantry rotations was done. Numerical experiments proved the following 
algorithms worked and demonstrated good image quality and high computational 
efficiency. Ramping down the speed of the moving table would inherently produce 

25 similar results. 
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The variable pitch (variable speed) spiral path C of the x-ray source is 
described by the following equations and depicted in Fig. 3, which shows 
mathematical notations of the spiral scan about the object being scanned: 

y^^Rcosis), jt 2 (j) = Asiii(j), = (1) 

5 Here 

s is a real parameter; 

z(s) is a function describing the third coordinate of the x-ray source on the 
spiral; the pitch is variable if z\s) is not a constant; 

R is distance from the x-ray source to the isocenter. 
10 The object being scanned is located inside an imaginary cylinder U of radius r , 
r < R (see Fig. 3). Let y/ be a smooth function with the properties 

^(0) = 0; 0<y/(0<l, /e[0,2*} (2) 
Even though it is not necessary, we will assume in addition that 

y/(0) = 0.5; ^ (2 * +l) (0) = 0, k>\. (3) 
15 Here and in what follows we assume that s Qi s l , and s 2 are always related by 

J|=^(^i-*o) +J o if s 0 <s 2 <s 0 +2n 9 (4) 

Ji^VOb-Jj+Ja if s 0 -27i<s 2 <s 0 . (5) 
Conditions (2) and (3) can be easily satisfied. One can take, for example, y/(i) = tl2 , 
and this leads to 

20 s x =(s q +s 2 )I2, s 0 -2n<s 2 <s 0 + 2n. (6) 

Denote 

*(W,)= ^-**» x t**>-**» sgn(, 2 -, 0 ) if 0<K-, 0 |<2*, (7) 
u(s KjoKW jf j — j m 
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Here 



y(s 0 )>y(S\)>y(s 2 ) are three points on the spiral related according to (4), (5); 



tt(j 0 ,^) is a unit vector perpendicular to the plane containing the points 



5 



y(s):=d 2 y/ds 2 . 

Any point strictly inside the spiral belongs to a PI segment A PI segment is a segment 
of line endpoints of which are located on the spiral and separated by less than one turn 
(see Fig. 4). We will assume that such a PI segment is unique. This holds, for 

1 0 example, if z \s) = const or if z \s) = const and z \s) does not change sign or if 
z \s) + z %y) does not change sign. Let s = s b (x) and s = s, (x) denote values of the 
parameter corresponding to the endpoints of the PI segment containing a 
reconstruction point x . We will call I FJ (x) := [s b (x),s t (x)] the PI parametric interval. 
The part of the spiral corresponding to I PJ (x) will be denoted C PJ (x) (see Fig. 4 

15 which illustrates a PI segment of an individual image reconstruction point). 
Next we fix a reconstruction point x inside the spiral and s 0 e I„ (x) . Find 
s 2 el pj (x) such that the plane through y(s 0 ),y(s 2 ) 9 and jK^Cv^)) contains x. 
More precisely, we have to solve for $ 2 the following equation 



20 Such s 2 exists, is unique, and depends smoothly on s 0 . Therefore, this construction 

defines s 2 := s 2 (s 09 x) and, consequently, u(s 0 ,x) := u(s 0> s 2 (s 0> x)) . Equation (9) can 

be solved by a variety of methods, all known under the name "root finders". The main 
reconstruction formula now is as follows: 



(*-J>C%))'*faA) = 0, s 2 eI FI (x). 



(9) 
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2n I x - y(s) | * oq I sin y 

(10) 

where 

/ is the function representing the distribution of the x-ray attenuation 
5 coefficient inside the object being scanned, 

e(s,x) = P(s,x)xu(s 9 x), 
x is the cross-product of two vectors, 
®(s,x,r) := cosrfi(s,x) +sinye(s,x) , 
D f is the cone beam transform of / : 

10 D f (y 9 &)=[ f(y + ®t)dt 9 

(11) 

P{s 9 x) := -j^fjf is the unit vector from the focal point y(s) pointing towards 
the reconstruction point x . 

1 5 Now we describe an efficient (that is, of the FBP type) implementation of inversion 
formula (10). It is clear from (9) that s 2 (s, x) actually depends only on s and /3(s,x) . 
Therefore, we can write 

u(s, P) := u(s y s 2 (s, P)) 9 e(s y £):=/? x u(s, 0), 0eS\ 

¥(*,£):= f ^D f (y(q%cmrfi+smre(sM ^—dy y 
* d 9 j^sin^ 

20 (12) 

/<*) = - A C (x U -, V(s,P(s,x))ds. 
2n ) inW\x-y{s)\ 

(13) 

11 
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Here S 2 is the unit sphere. 

To better understand equations (12), (13), we illustrate various important 
features on the detector array. Suppose first that the x-ray source is fixed at y(s) for 
some s e I PJ (x) . Project stereographically the upper and lower turns of the spiral 
5 onto the detector plane as shown in Fig. 5 which illustrates a stereographic projection 
from the current source position on to the detector plane used in the algorithm for the 
invention. 

Since the detector array rotates together with the source, the detector plane 
depends on s and is denoted DP( s) . It is assumed that DP(s) is parallel to the axis 

10 of the spiral and is tangent to the cylinder y 2 + = J? 2 (cf. equation (1)) at the point 

opposite to the source. Thus, the distance between y(s) and the detector plane is 2R . 

Introduce coordinates in the detector plane as follows. Let the d x -axis be 

perpendicular to the axis of the spiral, and the d 2 -axis - parallel to it. This gives the 

following parametric curves: 



l-cos(q-s) l-cos(q-s) 



h.<q-s<2n-L or &-2n<q-s<>-L, 
(14) 

where A is determined by the radius r of the imaginary cylinder U inside which the 
patient is located (see Fig. 3): A = 2cos\rlR) . The top and bottom curves are 
20 denoted and , respectively (see Fig. 6 which illustrates various lines and 
curves, such as boundaries, on the detector plane). The common asymptote of r tap 
and is denoted . Let x denote the projection of x . Since s e I„(x) , x is 
projected into the area between I\ and (see Fig. 6). Fix 
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s 2 e[s-27r + A 9 s+2x-Als 2 *s,mdlet U(s 2 ) denote the plane through 

7(^2) > ^ ?(*i(*,J 2 ))- If ^=^> n(s 2 ) is determined by continuity and 
coincides with the plane through y(s) and parallel to y(s), y(s) . The family of lines 
L(s 2 ) obtained by intersecting n(.s 2 ) with the detector plane is shown in Fig. 7. 
5 The main assumption under which equation (10) holds is that the curves T 

top 

and r fe , are convex. This happens, for example, if z\s) = const or if 

z \s) = const and z \s) does not change sign or if z '(s) + z m (s) does not change sign. 

By construction, given any x e U with 0(s,x) parallel to Il(s 2 ) and such that 
x appears to the left (rigfct) of the point of where the line L(s 2 ) intersects (T^) 
10 for the first time if x is above (below) 1^ , s 2 used here is precisely the same as ^ 
foimd by solving (9). The condition that we have formulated regarding the location of 
x relative to s 2 and 1^ guarantees that s 2 el pi (x). Since e(s,0)- 0 = Q 9 \e(s,P)\=l, 
we can write: 

fi = (cos0,sin0); e{s y 0) = (-sin0,cos0); frefrft e II(s 2 ). 
15 (15) 
Therefore, 




(16) 

Equation (16) is of convolution type and one application of Fast Fourier Transform 
20 (FFT) gives values of ¥(*, P) for all £ e n(^ 2 ) at once. 

Equations (13) and (16) would represent that the resulting algorithm is of the FBP 
type. This means that processing of every CB projection consists of two steps. First, 
shift-invariant and x -independent filtering along a family of lines on the detector is 

13 
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performed Second, the result is back-projected to update the image matrix. The 
main property of the back-projection step is that for any point x on the detector the 
value obtained by filtering at x is used for all points x on the line segment connecting 
the current source position y(s) with x . Since dldq in (16) is a local operation, each 
5 CB projection is stored in memory as soon as it has been acquired for a short period of 
time for computing this derivative at a few nearby points and is never used later. 
Now we describe the algorithm in detail following the six steps 10-60 shown in Fig. 2. 

Step 10. Load the current CB(cone beam) projection into computer memory. 
10 Suppose that the mid point of the CB projections currently stored in memory is y(s 0 ) . 
The detector plane corresponding to the x-ray source located at y(s 0 ) is denoted 
DP(s 0 ). 

Step 20. Fig. 8 is a four substep flow chart for identifying the set of lines, which 
corresponds to step 20 of Fig. 2. Referring to Fig. 8, the set of lines can be selected by 
15 the following substeps 21, 22, 23 and 24. 

Step 21. Choose a discrete set of values of the parameter s 2 inside the interval 

[s 0 - In + A, s 0 + In - A] . 

Step 22. For each selected s 2 compute the vector u(s Q ,s 2 ) according to 
equations (7), (8). 

20 Step 23. For each u(s 09 s 2 ) computed in Step 22 find a line which is obtained 

by 

intersecting the plane through y(s Q ) and perpendicular to the said vector 
u ( s o> s z) w tth * e detector plane DP(s 0 ) . 
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Step 24. The collection of lines constructed in Step 23 is the required set of 
lines (see Fig. 7 which illustrates a family of lines used in the algorithm of the 
invention). 

5 Step 30. Preparation for filtering 

Fig. 9 is a seven substep flow chart for preparation for filtering, which 
corresponds to step 30 of Fig. 2, which will now be described. 

Step 31. Fix a line L(s 2 ) from the said set of lines obtained in Step 20. 
Step 32. Parameterize points on the said line by polar angle y in the plane 
10 through y(s 0 ) and L(s 2 ). 

Step 33. Choose a discrete set of equidistant values y } that will be used later 
for discrete filtering in Step 40. 

Step 34. For each y } find the unit vector p } which points from y(s 0 ) towards 
the point on L(s 2 ) that corresponds to y s . 
1 5 Step 35. Using the CB projection data D f (y(q\ 0) for a few values of q 

close to s 0 find numerically the derivative {dldq)D f {y{q) 9 &)\ q ^ for all 

Step 36. Store the computed values of the derivative in computer memory. 
Step 37. Repeat Steps 31-36 for all lines L(s 2 ) identified in Step 20. This way 
20 we will create the processed CB data Wfo, ) corresponding to the x-ray 

source located at y(s 0 ) . 

Step 40. Filtering 
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Fig. 10 is a seven substep flow chart for filtering, which corresponds to step 40 
of Fig. 2, which will now be described. 

Step 41. Fix a line from the said family of lines identified in Step 20. 
Step 42. Compute FFT of the values of the said processed CB data computed 
5 in Step 30 along the said line. 

Step 43. Compute FFT of the filter 1/sin^ 

Step 44. Multiply FFT of the filter 1/siny (the result of Steps 43) and FFT of 
the values of the said processed CB data (the result of Steps 42). 
Step 45. Take the inverse FFT of the result of Step 44. 
1 0 Step 46. Store the result of Step 45 in computer memory. 

Step 47. Repeat Steps 41-46 for all lines in the said family of lines. This will 
give the filtered CB data <D(* 0 , /?,) . 

By itself the filtering step is well known in the field and can be implemented, 
for example, as shown and described in U.S. Patent 5,881,123 to Tarn, which 
15 is incorporated by reference. 

Step 50> Back-projection 

Fig. 1 1 is an eight substep flow chart for backprojection, which corresponds to 
step 50 of Fig. 2, which will now be described. 
20 Step 5 1 . Fix a reconstruction point x , which represents a point inside the 

patient where it is required to reconstruct the image. 
Step 52. If s 0 belongs to I PJ (x) , then the said filtered CB data affects the 

image at x and one performs Steps 53-58. If s 0 is not inside the interval 
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I PI (x) , then the said filtered CB data is not used for image reconstruction at 
x . In this case go back to Step 5 1 and choose another reconstruction point. 
Step 53. Find the projection x of x onto the detector plane DP(s 0 ) and the 
unit vector 0(s O9 x) , which points from y(s 0 ) towards x . 
5 Step 54. Using equation (9) identify the lines from the said family of lines and 

points on the said lines that are close to the said projection x . This will give a 
few values of <3>(s 09 /3j) for f3 } close to 0(s o ,x). 

Step 55. With interpolation estimate the value of <S>(s 09 fJ(s Q ,x)) from the said 
values of <b{s 09 fij) for fi s close to fi(s 09 x). 
10 Step 56. Compute the contribution from the said filtered CB data to the image 

being reconstructed at the point x by dividing <&(s O9 0(s O9 jc)) by 
-2**|*-j<Jk)l- 

Step 57. Add the said contribution to the image being reconstructed at the 
point x according to a pre-selected scheme (for example, the Trapezoidal 
15 scheme) for approximate evaluation of the integral in equation (15). 

Step 58. Go to Step 5 1 and choose a different reconstruction point x . 



Step 60. Go to Step 10 (Fig. 2) and load the next CB projection into computer 
memory. The image can be displayed at all reconstruction points x for which the 
20 image reconstruction process has been completed (that is, all the subsequent CB 

projections are not needed for reconstructing the image at those points). Discard from 
the computer memory all the CB projections that are not needed for image 
reconstruction at points where the image reconstruction process has not completed. 
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The algorithm concludes when the scan is finished or the image reconstruction 
process has completed at all the required points. 

The invention can work with other types of variable pitch(nonconstant speed) 
5 spiral scans. Fig. 12 shows an arrangement 500 of scanning an object 515 such as a 
human body, on a stationary table 510 within a spiral coil stand the object 5 1 5 being 
scanned remains stationary inside. The coil stand can be located inside of a chamber, 
or be a virtual coil stand within a chamber. As previously described, unlike the prior 
art, the invention is not limited to moving an object at a constant speed through a 

1 0 spiral scan. The object 5 1 5 can remain stationary within a stationary spiral coil type 
stand, where multiple x-ray sources SI, S2, S3, S4, S5, S6 and oppositely located 
detectors Dl, D2, D3, D4, D5, D6 arranged along the stationary coil stand 600 emit x- 
rays in a sequential manner about the stationary object 5 1 5 such as from right to left, 
left to right, the middle to the left, the middle to the right, and combinations thereof, 

1 5 to generate a spiral scan 

Still furthermore, the coil stand 600 can have fixed multiple x-ray sources and 
detectors so that the entire coil stand 600 can rotate about the object 515, and generate 
a spiral scan. 

Still furthermore, the spiral coil stand 600 can contain a single x-ray source SI 
20 and oppositely located detector Dl which moves along a spiral track on the stand 600 
about the fixed object 510 at constant and nonconstant speeds. Still furthermore, the 
spiral stand 600 can include coils links 610, 620, 630, 640, 650, 660, 670 that are not 
evenly spaced from one another so that the single x-ray source S 1 and opposite 
located detector Dl moving at a constant speed ends up passing along the length of 
25 the object 5 1 5 at different speeds. Thus, closely located links 6 1 0, 620 allow the 
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single source SI and detector Dl to pass at a slower rate over an object than distantly 
spaced apart coil links 650, 660, 670. 

The spiral coil stand embodiments described above can also work with 
constant pitch(constant speed) applications. 

5 

Other Embodiments of the invention are possible. For example, one can 
integrate by parts in equation (10) as described in the inventor's previous U.S. Patent 
Application 

Serial No. 10/143,160 filed May 10, 2002 now U.S. Patent 6,574,299, now 
10 incorporated by reference, to get an exact FBP-type inversion formula which requires 
keeping only one CB projection in computer memory. The algorithmic 
implementation of this alternative embodiment can be similar to and include the 
algorithmic implementation of Embodiment Two in the inventor's previous U.S. 
Patent Application 

15 Serial No. 10/143,160 filed May 10, 2002 now U.S. Patent 6,574,299, now 
incorporated by reference. 

Although the preferred embodiments describe applications of using x-ray 
sources for creating data for image reconstruction, the invention can be applicable 
with other sources such as but not limited to early arriving photons that create line 

20 integral data for image reconstruction. 

While the invention has been described, disclosed, illustrated and shown in various 
terms of certain embodiments or modifications which it has presumed in practice, the 
scope of the invention is not intended to be, nor should it be deemed to be, limited 
thereby and such other modifications or embodiments as may be suggested by the 

25 teachings herein are particularly reserved especially as they fall within the breadth and 
scope of the claims here appended. 
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I claim: 

1 . A method of reconstructing images from data provided by at least one 
detector, comprising the steps of: 
5 variable pitch spiral scanning an object with at least one detector with at least 

one cone beam projection source; and 

reconstructing an exact image of the scanned object in an efficient manner 
with a convolution based FBP (Filtered Back Projection) algorithm. 

10 2. The method of claim 1, wherein the scanning step further includes the step of: 
moving a table supporting the object through a rotational scanner with non-constant 
speed. 

3. The method of claim 1, wherein the step of reconstructing further includes the 
15 step of: 

shift invariant filtering of the cone beam projections; and 
back projection updating the image of die scanned object. 

4. The method of claim 1, wherein the step of reconstructing includes the steps 
20 of: 

storing approximately 2 to approximately 4 cone beam (CB) projections in 
memory at a time; and 

using one family of lines for the step of reconstructing. 

25 5 . The method of claim 1 , wherein the step of reconstructing includes the steps 
of: 
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storing 1 cone beam (CB) projection in memory at a time; and 
using one family of lines for the step of reconstructing. 



6. The method of claim 1, wherein the scanning step further includes the step of: 
5 supporting the object in a stationary position within a coil stand; and 

moving the one source and the one detector along the coil stand about the 
object at a constant speed to generate a spiral scan. 

7. Hie method of claim 1, wherein the scanning step further includes the step of: 
1 0 supporting the object in a stationary position within a coil stand; and 

moving the one source and the one detector along the coil stand about the 
object at a nonconstant speed to generate a spiral scan. 

8. The method of claim 1 , wherein the scanning step further includes the step oft 
15 supporting the object in a stationary position within a coil stand; and 

placing plural sources and oppositely located detectors on the coil stand about 
the object; and 

sequentially activating the cone beam projections to generate a spiral scan. 

20 9. The method of claim 1, wherein the scanning step further includes the step of: 
supporting the object in a stationary position within a coil stand; and 
placing plural sources and oppositely located detectors on the coil stand about 
the object; and 

rotating the coil stand about the object to generate a spiral scan. 

25 
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1 0. A method of computing images derived from a variable pitch spiral computer 
tomography with detectors, comprising the steps of: 

(a) collecting cone beam data from a detector during a variable pitch spiral 
scan of an object; 

5 (b) identifying lines on a planell intersecting the cone beam, wherein the step 

(b) of identifying lines includes the steps of: 

(bi) choose a discrete set of values of the parameter ,y 2 inside an interval 
containing s 0 , where .s 0 and s 2 are values of the parameter along the spiral; 
(bii) compute the vector u(s 0 , s 2 ) for each selected ^ 2 according to equations 
10 u(s iO- ^i)-^o))^fe)~yM) , n . , 0 



where, jis^yis^y^) points on the spiral related according to the 

1 5 following equations 

* = V(*2 " *o ) + *o if *o ^ *2 < % + 2^-, 
*i = VC*o -* 2 )+*2 if *o "2jt < tf 2 < 

and y/ is a function with the properties 
K0) = 0; 0<^(/)<l, /e[0,2*]; 
20 "(^o*^) i s a ^it vector perpendicular to the plane containing the 

points y(s 0 ly( Sl ) 9 y(s 2 ); 
y(s) :- dylds ; 
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(biii) find a line for each u(s Q9 s 2 ) which is obtained by intersecting the plane 

through y(s 0 ) and perpendicular to the said vector u(s Qy s 2 ) with the plane II; and 

(biv) repeating steps (bi-biii), and forming a family of lines from a collection 
5 of the lines ; 

(c) preprocessing and shift invariant filtering said data along said lines ; 

(d) back projecting said filtered data to form a precursor of said image; and 

(e) repeating steps a, b, c, and d until an image of the object is completed. 

10 11. The method of claim 10, wherein the scan includes an x-ray exposure of the 
object. 

12. Hie method of claim 10, further including the step of 

(b5) preparing the data prior to the step (c) filtering step of: 
15 (b5i) fixing a line L(s 2 ) from the set of lines obtained in step (b); 

(b5ii) parameterizing points on the lines by polar angle y in the plane through 
y(s 0 ) and L(s 2 ); 

(b5iii) choosing a discrete set of equidistant values y } \ 
(b5iv) find the unit vector P } for each y } which points from y(s 0 ) towards 
20 the point on L(s 2 ) that corresponds to y } ; 

(b5v) computer derivatives using the CB projection data D f {y{q) y &) for a 
few values of q close to s 0 using equation (d/dq)D f (y(q) y ®)\ q=So forail 
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(b5vi) store the computed derivatives; and 

(b5vii) repeat steps (b5i) to b5vi for all lines L(s 2 ) identified step (b), in order 
to create processed CB data "¥(s oy J3j) corresponding to the x-ray source 
located at y(s 0 ) . 

5 

13. Hie method of claim 10, wherein the back-projection step (d) includes the 
steps of: 

(di) fix a reconstruction point x > which represents a point inside the object 
being scanned where it is required to reconstruct the image; 
10 (dii) if s 0 belongs to I PI (x) , then the said filtered CB data affects the image at 

xand one performs Steps (diii) to (dviii), if s 0 is not inside interval I PI {x) , 
then the said filtered CB data is not used for image reconstruction at x and go 
back to step (di) and choose another reconstruction point Here I Pi {x) the PI 
parametric interval; 

15 (diii) find the projection x of x onto a detector plane DP(s 0 ) and unit vector 

J3(s 09 x) , which points from y(s 0 ) towards x ; 

(div) identify lines from family of lines and points on the said lines that are 
close to the said projection x , using equation 
(x-y(s 0 )) u(s 0> s 2 )^0 9 s 2 el„(x); 
20 (dv) using interpolate find value of O(s 0> 0(s o , x)) from ®(s 0 , fy) for pj 

close to fi(s 09 x) ; 

(dvi) compute contribution from filtered CB data to the image being 
reconstructed at the point x by dividing <S>(s 09 P(s oy x)) by -2x 2 \x-y(s 0 )\; 
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(dvii) add the contribution to the image being reconstructed at the point x 
according to a pre-selected scheme; and 

(dviii) go to step (di) and choose a different reconstruction point x . 

5 14. The method of claim 1 0, further comprising the steps of: 

storing approximately 2 to approximately 4 cone beam (CB) projections in 
memory at a time; and 

using one family of lines for the step of reconstructing. 

10 15. Hie method of claim 1 0, further comprising the steps of: 

storing 1 cone beam(CB) projection in memory at a time; and 
using one family of lines for reconstructing the exact image. 

16. The method of claim 10, further comprising the step of: 

1 5 moving the object at a nonconstant speed through a rotating source and 

oppositely located detector to generate the variable pitch spiral scan . 

1 7. The method of claim 1 0, further comprising the step of: 
positioning the object in a stationary position within a coil stand; and 

20 generating the spiral scan from at least one source and at least one oppositely 

located detector located on the coil stand. 

18. A method of identifying lines on a plane n used for reconstructing images 
based on a variable pitch spiral scan of an object in a computer tomography system, 

25 comprising the steps of: 
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(i) choose a discrete set of values of the parameter s 2 inside an interval 
containing s 0 , where s 0 and s 2 are values of the parameter along the spiral 
scan path ; 

(ii) compute the vector u(s 0 , s 2 ) for each selected s 2 according to equations 

3 w^ 0 ,^ 2 ;- — — — sgn(j 2 -^ 0 ) if 0< hy 2 ~^ 0 <2;r, 

IW^)-J<^))xO<^)-J<Jb))l '2 01 



where, yte^yfalyisj are three points on the spiral related according to the 
1 0 following equations 

Si=V(s 2 -So)+s<> if * 0 ^2<*o+2;r, 
Si=Hs Q -s 2 )+s 2 if s 0 -2x<s 2 <s 0 , 

and y/ is a function with the properties 
y<0) = 0; 0<v/(0<l, / e [0,2*] ;w(.y 0 ^ 2 ) is a unit vector 
15 perpendicular to the plane containing the points y(s Q \y(s x \y(s 2 ) ; 

y(s)v=dylds\ 
y(s):=d 2 y/ds 2 ; 

(iii) find a line for each u(s 0 ,s 2 ) which is obtained by intersecting the plane 
through y(s 0 ) and perpendicular to the said vector u(s 0 ,s 2 ) with the plane II ; and 
20 (iv) repeating steps (i-iii), and forming a family of lines from a collection of 

the lines so that an image is reconstructed using die lines. 
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19. Hie method of claim 18, further comprising the step of: 

moving the object at a nonconstant speed through a rotating source and 
oppositely located detector to generate the variable pitch spiral scan . 

5 20. The method of claim 1 8, further comprising the step of: 

positioning the object in a stationary position within a coil stand; and 
generating the spiral scan from at least one source and at least one oppositely 
located detector located on the coil stand. 
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Fig.2 



Step 10. Load the current CB (cone beam) projecticn into computer 
memory. Suppose that the mid point of the CB projections currently 
stored in memory is y(s 0 ) . 



± 

Step 20, Identify the fenrily of filtering lines. Then, choose a discrete set 

of lines from that family. 



Step 30, Preparing for filtering. Parameterize points on the said lines 
selected in Step 30 by polar angle/ . Using interpolation, compute the 
derivative of the CB data (dfdq)D f (y(q) 9 0) at points on the said 
lines that correspond to discrete values of the polar angle. 



1 


r 


Step 40. Filtering. For each line identified in Step 20 convolve the data 
for that line computed in Step 30 with filter 1/siny . 


J 


I 



Step 50- Back-projection. For each reconstruction point x, back- 
project the filtered data found in Step 40 according to equation (10). 
Then goto Step 10, unless there are no new CB projections to process 
or image reconstruction at all the required points x have been 
completed. 
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Fig. 5 
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Fig.8 



Step 20. Finding femilies of lines for filtering 



Step 21. Choose a discrete set of values of the parameter s 2 inside 


the interval fo-^+A^ +2/r- 


-A]. 




t 


Step 22. For each selected s 2 compute the vector u(s 0 .s 2 ) 


according to equations (7), (8). 






f 



Step 23. For each u(s 0 ,s 2 ) computed in Step 22 find a line which 
is obtained by intersecting the plane through y(s 0 ) and 
perpendicular to the said vector u(s 0 , s 2 ) with the detector plane 
DP(s 0 ). 



Step 24. The collection of lines constructed in Step 23 is the 
required set of lines (see Fig. 7 which illustrates a family of lines 
used in the algorithm of the invention). 

Step 30 
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Fig. 9 

Step 30. Preparation for 
filtering 

Step 31. Fix a line L(s 2 ) from the said set of lines obtained in 
Step 20. 



Stej>32. Parameterize points on the said line by polar angle y in 
flie plane through y(s 0 ) and L(s 2 ) . 

— — 2: 

Step 33. Choose a discrete set of equidistant values y } that will 

be used later for discrete filtering in Step 40. 



Step 34. For each y J find the unit vector fij which points from 
y(s 0 ) towards the point on L(s 2 ) that corresponds to y s - 



Step 35. Using the CB projection data D f {y(q),&) for a few 
values of q close to s 0 find numerically the derivative 
(d/dq)D f <j(q) 9 e)l^f6rall & = 



Step 36. Store the computed values of the derivative in computer 
memory. 



Step 37. Repeat Steps 31-36 for all lines L(s 2 ) identified in Step 
20. This way we will create fee processed OB data ^(s^^) 

Step 40 
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Step 41. Fix a line from the said family of lines identified in Step 
20. 



Step 42. Compote EFT of die values of the said processed CB data 
computed in Step 30 along the said line. 



Step 43. Compute FFT of the filter 1 /siny 

i 

Step 44. Multiply FFT of the filter 1/siny (the result of Step 43) 
and FFT of die values of the said processed CB data (the result of 
Step 42). 

i 

Step 45. Take the inverse FFT of the result of Step 44. 



Step 46. Store the result of Step 45 in computer memory. 



Step 47. Repeat Steps 41-46 for all lines in the said family of 
lines. This will give the filtered CB data fl>(s 0 >/?,) . 

I 

Step SO 
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Fig. 11 

Step 50* Back-projection 

Step SI . Fix a reconstruction point x , which represents a point inside 
the patient where it is required to reconstruct the image. 

I 

Step 52. If s 0 e I„(x), then the said filtered CB data affects die image 
at x and one performs Steps 53-58. If s 0 € I„(x) , then the said filtered 
CB data is not used for image reconstruction at x . In this case go back 
to Step 51 and choose another reconstruction point 



i 


Step 53. Find the projection jc of x on 


to die detector plane DP(s 0 ) 


and the unit vector fi(s oy x) , which poi 


ntsfrom y(s 0 ) towards x. 



I 



Step 54. Using equation (9) identify the lines from die said family of 
lines and points on the said lines that are close to die said projection x . 

i ■ 

Step 55. With interpolation estimate the value of O(^ 0> /9(j 0 ,x)) ftom 
the said values of <S>(s„Pj) for p } close to fi(s 0> x). 

I 

Step 56. Compute the contribution from the said filtered CB data to the 
image being reconstructed at the point x by dividing <P(s 0 ,/3(s 0 , x)) 

fry-2»*|x-jfr 0 )l. 

i 

Step 57. Add the said contribution to die image being reconstructed at 
the point x according to a pie-selected scheme for approximate 

evaluation of the integral in equation (15). 

v 

Step 58. Go to Step 51 and choose a different reconstruction point jr. 
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